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Abstract 

We study the influence of intracluster large scale magnetic fields on the thermal 
Sunyaev-Zel'dovich (SZ) effect. In a macroscopic approach we complete the hydro- 
static equilibrium equation with the magnetic field pressure component. Comparing 
the resulting mass distribution with a standard one, we derive a new electron density 
profile. For a spherically symmetric cluster model, this new profile can be written 
as the product of a standard profile and a radius dependent function, close to 
unity, which takes into account the magnetic field strength. For non-cooling flow 
clusters we find that the observed magnetic field values can reduce the SZ signal 
by ~ 10% with respect to the value estimated from X-ray observations and the 
/3-model. If a cluster harbours a cooling fiow, magnetic fields tend to weaken the 
cooling flow influence on the SZ-effect. 
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1 Introduction 



The SZ-effect is rapidly turning into an important astrophysical tool thanks to 
the progress of the observational techniques, which allow increasingly precise 
measurements. In view of these developments it is thus relevant to study fur- 
ther corrections to it, such as relativistic effects ( |Kephaefi| , |1995D , the shape of 
the galaxy cluster and its finite extension or a polytropic temperature profile 
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(see e.g. |Puy et aT] (pOOOj )), corrections induced by halo rotation ( |Cooray and 
Chen| , |2001|) , Brillouin scattering ([Sandoval- Villalbazo and Maarten^ , 2002 ), 



early galactic winds ( [Majumdar et al.] , |2001| ) and the presence of cooling flows 
(|Schlickeiser| , |1991| ; [Majumdar et al.| , [200 1[) . These additional effects are of 
different relevance and often depend on the specific cluster values. 



Whereas, e.g. cooling flows are not present in every cluster of galaxies, the need 
for relativistic corrections due to energetic non-thermal electron populations 
seems to be common in most clusters. ( [Blasi et "al] , [2000[ ; [Shimon and Rephaeli| , 
2002| ). These relativistic electrons produce a hard X-ray component in excess 
of the thermal spectrum by Compton scattering off the Cosmic Microwave 
Background (CMB) and by non-thermal bremsstrahlung. Their emission has 
been quite possibly detected in Coma ( Kephaeli et al. 



1999[; jbusco-i'emianc 



"eTaTl , [19991) , A2199 ( [Kaastra et ^ [1999D , A2256 ([Molendi et al] , |2000D and 



A2319 ([Gruber and Kephaeli[ , pOO^ ) by RXTE and BeppoSAX satellites. The 
main evidence for the existence of relativistic electrons comes from radio syn- 
chrotron emission of extended intracluster regions ( [Giovannini et aT] , [1999| ; 



Giovannini and l^eretti , 2000| ) and is, therefore, closely related to the presence 
of magnetic fields. The magnetic fields in the intracluster gas lead to accel- 
eration processes and modify the classical Maxwell-Boltzmann distribution of 
the electrons, which might acquire a significantly non-thermal spectrum and 
thus account for the observed hard X-ray spectra ( [Ensslin et al. 



T999[; [BTiiij 



2000| ) . Consequently, several authors ( |Rephaeli| , [1995[ ; [Itoh et al\ , [1998[ ; [Chall: 



nor and Lasenby[ , [1998[ ; [Birkinshaw[ , [1999[ ) derived relativistic corrections to 



the thermal SZ-effect up to different leading orders. Though the magnetic field 
is related to the relativistic electron population, its own influence - as a non- 
thermal cluster component - on the SZ-effect has not yet been investigated. 
Whereas the importance of relativistic corrections to the SZ-effect depends on 
the cluster temperature, the magnetic field seems to be ubiquitous with a mean 
field value of 5 — 10 fiG in the cluster cores ([Clarke et al.[ , [200 1[) . Independently 
of the diffuse non-thermal radio emission, excess Faraday rotation measure of 
polarized radio emission in radio sources within or behind the cluster can 
prove the existence of magnetic fields. This method was applied to the Coma 
cluster, where [Feretti et al.[ ([1995[) found a large magnetic field oi B > 8.3 fiG. 
[Govoni et al. (1999) found a similar value (5 — 10 fiG) for A119. Clarke et al 



( [1999[) derived magnetic field strengths of a few iJ,G for their cluster sample 
by using a statistical Faraday rotation measure technique. Unfortunately, the 
structure of the magnetic field is presently poorly known. Contrary to the 
static situation in non-cooling flow clusters, the magnetic field is believed to 
become dynamically significant in the cores of cooling flow clusters ( Eilek an^ 
Owen[ , |2001|) . [Taylor et al.j ([2001[) found a magnetic field strength of up to 
40 fiG in the Centaurus cluster. The converging cooling flow (see e.g. [I'abian 



et al.[ ([1991[) ; [Fabian[ ( [1994[) ) causes a compression and enhancement in the 



magnetic field strength and finally reconnection might transfer the magnetic 
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energy back to the plasma when the magnetic field pressure becomes compa- 
rable to the thermal gas pressure. 

As the ratio of the magnetic pressure (Pg) to the gas pressure (Pg) reaches 
^ ~ 10~^ for the quoted values of non-coohng flow clusters and even unity 
for the cores of coohng flow clusters, we consider the additional magnetic held 
pressure term to be significant and we will examine its infiuence on the SZ- 
effect. We remark that all the quoted values refer to large scale magnetic fields 
with a coherence length of typically 1 — 10 kpc. There are essentially no useful 
limits on the strength of any small scale magnetic fields. Our calculation is 
based on a macroscopic picture inferred from the existence of this additional 
(large scale) pressure component and we do not start our considerations at the 
level of the single particle movement in a magnetic field. The current cluster 
data reveal magnetic field values which result to be significant for a correct 
SZ analysis, especially towards the cluster core where the electron density in- 
creases. 

The aim of the paper is to examine the infiuence of large scale magnetic fields 
on the thermal SZ-effect. We, therefore, choose a phenomenological approach, 
where the addition of a magnetic field pressure to the gas pressure is well 
justified. For a given magnetic field model, we can then estimate the change 
in the electron density and the temperature profiles as compared to standard 
ones used in the literature in the absence of magnetic fields. 

The paper is organised as follows: In section 2 we present the theoretical 
model for the magnetic field contribution. We derive new gas density profiles 
which are then used to calculate the SZ-effect. We distinguish two situations 
according to whether a cluster harbours a cooling fiow or not. Section 3 shows 
our results and contains a discussion of how magnetic fields influence the SZ 
signal and what are the observational consequences. As an illustration we 
apply our results to the non-cooling flow cluster A119. Our conclusions are 
given in section 4. 



2 Magnetic field contribution 

2.1 Non-cooling flow dusters 

We suppose a standard spherically symmetric model for a cluster, which is 
assumed to be in a relaxed state with a static gravitational potential. Thus, 
hydrostatic equilibrium can be expected. The intracluster plasma is treated 
as an ideal gas and thus the well known hydrostatic equilibrium equation has 
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the form (see e.g. ^arazin] (|1988| )): 

1 dPgjr) 
Pg{r) dr 



GMir) 



where G is the gravitational constant, Pg{r) and Pg{r) are the radius dependent 
gas density and pressure, respectively, and 7W(r) is the gravitating mass within 
the radius r. M.{r) is mainly determined by the dark matter profile, whereas 
the gas mass contribution is negligible. When taking into account magnetic 
fields, the hydrostatic equilibrium Eq.(l) has to be completed with a magnetic 
hydrostatic pressure term Pb of the form: 

Psir) = (2) 



where B{r) is the magnetic field strength, which is supposed to be spherically 
symmetric. Thus, the magnetic field contributes to the total pressure opposing 
the gravitational force. The gas pressure Pg{r) in Eq.(l) is replaced by the sum 
of gas and magnetic field pressure: 

P,(r) + Psir) = -^p,{r)T,{r) + (3) 



where k is the Boltzmann constant, p the mean molecular weight, the 
proton mass and Tg{r) the gas temperature at radius r. 

When deriving the mass distribution from Eq.(l), either with only the gas 
pressure or with the magnetic field pressure added as in Eq.(3), we obviously 
end up with the two different^ equations: 



M{r) 



G firrip 
kTsir) 
G finip 



dlnpg{r) d\nTg{r 

dlnr 
dlnpB{r) 



dlnr 



+ 



dlnr 
d In Tb (r 
dlnr 



dPBir 



G pb{t) dr 



(4) 
(5) 



where Tsir) and pB(r) describe the gas temperature and density, respectively, 
in the presence of the magnetic field. 



As next we compare the Eqs.(4) and (5). To find the cluster mass distribution 
we can either proceed our analysis according to Eq.(4) or, if we take into ac- 
count magnetic fields, according to Eq.(5). However, the true total gravitating 

^ If necessary we label quantities with magnetic fields with an index 'B' to distin- 
guish them from quantities without magnetic fields. 
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cluster mass is unique and must be the same, thus Aisiri) = A4{ri), ri being 
the cluster limiting extension. Nevertheless, Tsir) and p_B(r) can be different 
from Tg{r) and Pg{r) due to the magnetic field pressure gradient^. In what 
follows, we want to relate psir) to Pg{r), which is supposed to be known and 
to be a correct theoretical description for the gas density profile in absence of 
magnetic fields (e.g. /3-model, |Sarazin| (|1988|) ). 



As dark matter is the dominant mass component in clusters of galaxies (up 
to 80 — 90% of the total mass) and supposed not to be affected by magnetic 
fields, we expect in good approximation the equality between the Eqs.(4) and 
(5) to be satisfied for all radii r: 

MB{r)=M{r). (6) 



Clearly, this approximation is justified as long as the density profiles pB(r) and 
Pg{r) do not differ substantially. As we will see, our results do indeed satisfy 
this requirement. For the sake of simplicity we assume isothermal temperature 
profiles: Tg{r) = Tg, Tsir) = Tb- By setting equal the right hand sides of the 
Eqs.(4) and (5), we find a first order differential equation for pB{r), which 
yields: 



pB[r) 



Tb 



' Pg,0 



1 + P 



9,0 



-i^ pmp \ 



P'si 



Pgir) 



■ dr 



(7) 



where prime is the derivative with respect to r and pg^ is the central gas den- 
sity. The boundary condition is chosen to be psiji) = Pgiji) for r/ at the cluster 
limiting radius, where the magnetic field is negligible. Physically, we do not 
expect the temperatures Tg and Tb to differ significantly. Thus, in the sense of 
a first order development, we suppose Tg and Tb to be approximately equal, 
which simplifies the above equation. Mathematically, the boundary condition 
and the assumption of isothermal temperature profiles even impose Tb = Tg-. 
At the limiting radius r^, the magnetic field pressure gradient in Eq.(5) van- 
ishes and the densities are then equal. As the temperatures are isothermal, 
the second term in brackets in the Eqs. (4) and (5) drops, stating that Tg = Tg 
at ri, which is then also true for the whole cluster. 



To evaluate Eq.(7) we need a magnetic field model. Various physical models 
have been invoked to explain the rotation measure in radio sources in clus- 
ters of galaxies ( |J aff'e| , p.980| ; [Iribblel , |1991D . According to them, the magnetic 



^ We note that our starting point is different from the usuai observer's point of 
view: From a data set, a density profile is first derived and later the magnetic field 
contribution is added as in Eq.(5). This procedure gives a higher cluster mass (see 
e.g. iLoeb and Mao| (|1994D ). 
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field distribution is correlated with the electron density of the thermal gas. 
Recently, Polag et al.| ( 2001 ), based upon a correlation between X-ray surface 
brightness and Faraday rotation measure, derived the following relation: 



B{r) oc (ne(r)) ' 



where ne(r) is the electron number density and 7 the slope of the B — Ue 
relation, which depends on the specific cluster of galaxies. They showed, that 
Eq.(8) can be motivated by both simulations and observational data. Fur- 
thermore, they clearly excluded the possibility of a constant magnetic field 
through the intracluster medium. With the relation (8) the magnetic field 
profile is proportional to PB(r) and we can thus rewrite Eq.(7) as follows: 



pB{r) = pg{r) 



1 + 



1 5, 



27 
Pbo 



-( 



11 nip 
87r V kTa 



Pgir) 



dr 



(9) 



where pb,o and Bq are the cluster central gas density and the central magnetic 
field value, respectively. Eq.(9) is an integro-differential equation for psir). 
It expresses the modified density psi^) as the standard density pg{r) from 
Eq.(4), multiplied with a radius dependent function which involves the mag- 
netic field pressure. As in perturbation theory, pB(r) can also be interpreted 
as the product of the unperturbed quantity Pg{r) multiplied by the function 
in brackets in Eq.(9), which is close to unity. Following an iterative procedure 
to solve Eq.(9), Pg{r) is reinserted into the function in brackets. We stop after 
the first iteration to get: 

psir) ^ pg{r) 



1 + 



27 ■ 



Bl 



ml 



Pgij) 



dr 



(10) 



Setting Pg{r) = Pgflf{r), where /(r) is the shape of the gas profile, Eq.(lO) 
finally reads: 



Psir) =pg^of{r) 



9 "-gfi' 



fir) 



dr 



< Pg{r) 



:iii 



where the last inequality arises immediately because the magnetic field pres- 
sure decreases towards the cluster boundary. The modified gas density pb can 
thus be calculated from any standard density Pg, the cluster temperature Tg 
and some central magnetic field value Bq. We stress that this result follows 
from the starting point that the cluster mass A4{r) = Aisir) can be deter- 
mined in two different ways following the Eqs.(4) and (5), when the self-gravity 
of the gas is neglected. 
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2.2 Cooling flow clusters 



We suppose again spherical symmetry. For simplicity, we adopt a homogeneous 
steady-state cooling flow model. The gas has a single temperature and density 
at a given radius and no mass drops out of the flow. The cluster is expected 
to be in a relaxed state, so that hydrostatic equilibrium allows us to use an 
isothermal (3-mode\ (|Sarazin| , |1988|) . The dynamics in the cooling flow region 
can thus be described by a set of Euler equations. Mass, momentum and 
energy conservation read ( Mathews and Bregman , 1978 ; White and Sarazin , 
rW|,|; |Sarazin| , [1^ : 



^^('"VsW^lO) = 0, 

, ,dv(r) 1 dPJr) GM(r) 



dr 



Pg{r) dr 



dE{r) Pg{r) dpg{r) 



dr 



p2(r) 



dr 



0, 



-Apg{r). 



(12) 
(13) 

(14) 



Here r, v{r) and Pg{r) are the radius, gas mean velocity and gas pressure, 
respectively in the cooling flow. The velocity v{r) is deflned to be negative 
for the inward directed cooling flow. The internal energy is E{r) = {3/2)6{r) 
with the temperature parameter 6, which deflnes the square of the isothermal 
sound speed c^: 



e(r) := c^(r) = 

finip 



(15) 



where Tg is the gas temperature in the cooling flow. Ai{r) is the gravitating 
cluster mass inside the radius r. In order to determine it, the above mentioned 
authors used the following dark matter {DM) density proflle: 

pDM{r) = f° (16) 
1 + [r/rc) 



with a central cluster density po = 1-8 ■ 10~ g cm'' and a core radius = 
250 kpc. As usual, we assume that the cooling flow makes no signiflcant con- 
tribution to the cluster mass density and that the gas self-gravity can be 
neglected. 

The cooling function A(^) is deflned so that Ap^ is the cooling rate per unit 
volume in the gas. We use an analytical flt to the optically thin cooling func- 
tion (pR^aymond et al.| , |1976| ) as given by ( ^arazin and White] , |1987| ; [Majumdai 
and NathI , 120001) : 
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MO) 



IQ-'^'^erg cm? 

exp 



+ 0.313 -T^ "^ 



4.7 • exp 
T 



T 



3.0 ■ IQ^K 



3.5 ■ lO^fs: 

4.4' 



4.5' 



+ 6.42 ■ T-°-2 ■ exp 
+ 0.000439 -T^-^l 



T 



2.1 ■ lO^/sT 



4.0- 



(17) 



As noted by |Majumdar and Nath| (|2000| ), this fit is accurate to within 4% for 
a plasma with solar metalicity in the temperature range 10^ < T < 10* K. 
For 10* < T < 10^ K, it underestimates cooling by a factor of order unity, 
compared to the exact cooling function as in [Schmutzler and Tscharnuter 
(|1993|) . The continuity Eq.(12) directly yields an expression for the gas density 
Pg{r) in the cooling flow: 



m 



Pg[r) 



47rr2 v{r) 



where rh < is the constant cooling flow mass deposition rate which enters 
as a parameter in our model. The Eqs.(12)-(18) describe our standard cooling 
flow model without magnetic flelds. 



To include the magnetic flelds in our calculation we follow the paper by [Soker 



and Sarazin] ( |1990| ). As it was already argued by them, we also limit our 
discussion to large scale magnetic flelds. We assume that the magnetic fleld 
lines are frozen-in to the inward flowing homogeneous cooling gas. Under this 
assumption, the radial and tangential coherence lengths vary as (poker and 
Sarazin 



19901 ): 



cool 



' VB{r) ' 

. VBrOol . 



cool 



^cool 



(19) 



where Icooi and vb,cooI are the typical coherence length and inflow velocity, re- 
spectively, at the cooling radius Vcooi ~ 100 kpc. The magnetic fleld is assumed 
to be isotropic outside of rcooi, so that Ij. = k = hooi and = /2 = B^^^JS 
for r > Tcooh where Bcooi is the magnetic fleld strength at Tcooi- Inside the 
cooling flow region the radial and the tangential fleld components are then: 

BAr) - {^)\ BM ^ f-B,^. (If) (^) . (20) 
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The compression of the gas is expected to produce a sensible increase of the 
frozen-in magnetic field strength. Thus, the magnetic field lines become in- 
creasingly radial as the gas flows inward. As it is pointed out by pitti et aT 
(|2002|) , this corresponds to the physical condition in the cooling flow with 
vt <^ If^l, where vt and vb are the turbulent velocity and the mean inflow 
velocity, respectively. In this model, the turbulence does not disturb the field 
geometry during infall. The case vt 3> \vb\ is considered by |Tribble| ( [19931 ). 
The field strength in the Soker and Sarazin model grows faster and in the 
central region the field reaches values which are an order of magnitude larger 
than the one of Tribble, which becomes more isotropic towards the center. As 
it is mentioned by [Gitti et aLl (|2002|) , there seems to be observational evidence 
for Tribble's model in the Perseus cluster, where it was found vt ~ 60 km /s 
and [fsl ~ 20 km /s at Tcooi- Contrary to this result, different authors (|Clarke| 
et al.| , |1999| ; [I'aylor et al.| , |1999| , |2001| ) quote magnetic field values of the order 
B ^ 10 — 100 fiG for the cores of cooling flow clusters. Furthermore, intense 
magnetic fields of the order of some tens of fiG have been derived by detection 
of extremely high Faraday rotation measures throughout radio galaxies in the 
centers of cooling flow clusters ( |Ge and Owen| , |1993| ; [Taylor and Perley| , [1993[ ) . 
These results would favour the model by Soker and Sarazin. Being aware of 
this discrepancy, we will adopt the Soker and Sarazin model for the following 
presentation. The magnetic field pressure can then be expressed from Eq.(20): 



PB(r) 



B,cool 




(21) 



where Pb,cooI is the magnetic field pressure at the coohng radius Vcooi- From 
their discussion about magnetic forces and small scale magnetic fields, pokei 
and Sarazin[ ( [1990[ ) concluded in adding simply the magnetic field pressure 
term to the gas pressure in the Eqs.(13) and (14). [Gonqalves and Friaqa[ (|99|) 
used then the same approach for their simulation of the evolution of the intr- 
acluster medium with magnetic field pressure. 

Proceeding like this and eliminating the cooling flow gas density Pb{i^) from 
the Euler Eqs.(13) and (14) with Eq.(18), we end up with a system of two 
coupled first order ordinary differential equations for the isothermal sound 



speed Cs^b{^) = \/0b{i^) and the infall velocity ^^(r) in r: 



dVB 

dr 

dOB 
dr 



Vb ^GM - lOrO 



B 



2-K vi 



B,coolfl 



r\WB - 3vl + 8^PB,cooif: 



(22) 



eB{2rvl-GM) + 



m 
i-K vi. 



B 
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B 



l^-PB.cooi/a] + ^^PB,cool9Bf2 



{^Ob - 3vl + 8^PB,coolf2 
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m 6 \ 



These equations reduce to the system of equations derived by [Mathews and 



Bregmah] ( |1978| ) in the absence of magnetic fields, Pb,cooI = 0. The functions 



fi are defined as follows: 



3 r ?> Vb ' 

4 2 2 

f I \ '^B'^cool cool^B, cool f^r\ 

72 r = — , 25 

r Vb 

f (^\ — '^B^cool I 8 ^cool'^B,cool 

Sr-^ 3 Vb 

/^(^) = (^!i^y + 2(^^)'(^^^^)'. (27) 



Both Eqs.(22) and (23) have singularities at the sonic radius r^. Under favor- 
able conditions, the sonic singularity corresponds to a crossing of two critical 
solutions which are asymptotes to families of hyperbolae near the singular- 
ity. This was discussed by [Mathews and Bregman] ( [197^ ) in the absence of 
magnetic fields. The Eqs.(22) and (23) allow for transitions from subsonic to 
supersonic flows, if the numerators and denominators vanish at r^: 



= 59B{rs)-3vl{rs) + S-PB,cooif3{rs), (28) 

m 

( 9 , , 167r ^ p / x\ J^AfTo) 
= r.(64(r.)-— PW.W)+^^ 

+—PB,cooifi{rs) - 3GA<(r,), (29) 
m 

but the quotients are well behaved. Given the form for A^(r), A(T) and the 
parameters Pb,cooI, Vcooi and vb,cooI, we could solve these equations for possi- 
ble sonic point (s) r^, depending on the chosen temperature Tg. This procedure 
was adopted by several authors who described the cluster dark matter density 
with a King profile. As it was noted by [Mathews and Bregman[ ( [1978| ) and 
[Sulkanen et al.[ ( p.989| ), this profile allows either one or three solutions for in 
Eq.(29) in the absence of magnetic fields. The sonic radius is then given by 
the smallest value. To find solutions for the Eqs.(22) and (23), in the absence 
of magnetic fields, the above authors started the integration away from the 
sonic point towards the cooling radius Tcooi- Since the expressions for the 
derivatives are indeterminate at r^, they had to be replaced by nonsingular 
expressions, derived by making use of the BernouUi-de FHopital's rule, which 
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have to be matched to hydrostatic equihbrium outside the coohng flow region. 
As the presence of magnetic fields complicates substantially this procedure, 
we avoid this time-consuming method. For simplicity, we choose a different 
approach which we describe in section 3.2. The magnetic field influence in the 
region outside rcooi ~ 100 kpc is calculated by the method described in section 
2.1. 



3 Results and discussion 



3.1 Non- cooling flow clusters 



The scope of this section is twofolded: 
pact of the modified density psir), and 
consequences and apply as an example 
A119. Based on the new density profile 
SZ-effect. 



First, we discuss the theoretical im- 
second, we look at the observational 
our calculations to a specific cluster: 
Pb{j), we calculate the change in the 



3.1.1 Theoretical considerations 

King ( |1966| ) developed a self-consistent truncated density distribution and 
gave an analytic function for a cluster potential. Based on the hydrostatic 
equilibrium Eq.(l) and a constant temperature, one derives then the well 



known isothermal /?-model for the gas density (for a review, see e.g. ^arazin 



( |1988| )). Since then, this model has been extensively used in the literature. In 



the present discussion we will limit the standard input density Pg in Eq.(ll) 
to this /5-model: 

p,(/?,r) = p,,o/(r) = p,,o 1+ -) , (30) 



where is the cluster core radius and (3 a fitting parameter. As the /5-model is 
an appropriate theoretically motivated description in the absence of magnetic 
fields, Eq.(ll) states that, observational data, which involve magnetic fields, 
should not be fitted with a /3-model. The correct fitting function is related 
to a /3-model, Pg{P,r), through the function in brackets in Eq.(ll): 

PBir)=Pg{f3,r)[l + hB{r)], (31) 



where we introduced the notation /^^(r) for the radius dependent magnetic 
field contribution. Eq.(31) clearly shows, that the gas density profile p_B(r), as 
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Fig. 1. The modified profile psif) compared to the /3-profile pg{f3,r) for standard 
cluster values. The mean molecular weight is p = 0.63. The profiles are normalized 
by the central gas density pg^. 

used for the surface brightness Sx or energy spectra, is no longer of the type 
of a /3-model. 

To get an estimate on how significantly psir) varies compared to Pg{P, r), we 
adopt some standard cluster values, which should be free of the magnetic field 
influence. The use of these standard (observational) values will be justified in 
more detail in the next section 3.1.2. The input parameters are: Tc = 250 kpc, 
Tg = 2- 10^ ir, Ugfl = 1.2 ■ 10-2cm-3, (3 = 2/3, Bo = 7.5 pG and 7 = 0.9. For 
the isothermal /3-model Pg{P,r), as introduced in Eq.(30), the function hsir) 
can be found analytically to yield: 



hB(r) = A- 



27 
27 - 1 



1 + 



3/3(1/2-7) 



1 + 



3/3(1/2-7) • 



(32) 



with A = $- — ), 7 ^ 0.5. 

Fig. 1 shows the normalized density profiles for pg{(3,r) and psif)- We note 
that the modified profile Pb{t) is lower than Pg(/?, r) with the biggest difference 
of 10 — 20% in the most inner part of the cluster, where the magnetic field 
becomes important. 

Based on these profiles, the SZ-effect can be calculated. The frequency depen- 
dent intensity change of the cosmic microwave background photons (CMB) 
due to inverse Compton scattering off the hot intracluster electrons can be ex- 
pressed as follows ( [Sunyaev and Zerdovich] , [1972^ [Rephaeli] , |1995| ; [Birkinshaw , 
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19991 ): 



AJ(x) = iog{x 




dL 



(33) 



where x = |^ is the dimensionless frequency with T the CMB temperature 

and zq = ^|^^2 • The function g{x) defines the spectral shape of the ther- 
mal SZ-effect. The integral is the Comptonization parameter y describing the 
cluster properties with Tg,me the electron cluster temperature and electron 
mass, respectively. Ue is the electron number density in the cluster and ctt 
the Thomson cross section, {k, h, c are the Boltzmann constant, the Planck 
constant and the speed of light, respectively.) As we are interested in the in- 
fluence of magnetic fields, we define a to be the ratio between the SZ-effect 
with the modified density pB(r) and the standard SZ-effect with pg{P,r). For 
our standard values we find: 

a:=^ = -p^fPZ-- = 0-92, (34) 
A/ j;;' pg{(3,r)dr 



where the integration is along a line of sight through the cluster center with 
a limiting radius ri = 1000 kpc. From our theoretical analysis we thus find 
that magnetic fields reduce the SZ-effect by ~ 10% compared to the value 
one would get for an average standard cluster without magnetic fields. 
Indeed, a = a{Bo,Tg,ngfi,ri, jS,^) depends on several observational values. 
Modifying the central magnetic field strength Bq, while keeping constant the 
other parameters, can substantially change the ratio a. Bq ^ 1 pG reduces 
the SZ-effect by less than 1%, whereas Bq ^ 10 pG gives a = 0.85. Similarly, 
a lower temperature Tg = 10^ K combined with our standard values gives 
a = 0.84, whereas a very hot cluster with Tg = 10^ leads to a = 0.98. 
In an analogous way, a decreases with a smaller central gas number density 
Ugfl. Reducing rip q by 20% gives a = 0.90, instead a 20% increase results in 
a = 0.93. A larger cluster limiting radius ri makes the magnetic field influence 
generally more important, a is not sensitive to little changes in the parameter 
P, because both Pg{P, r) and pBi^) depend on it. Polag et "aL| ( p001|) estimated 



the slope 7 for the n^-B relation to be in the range 7 G (0.5 — 1.0). a decreases 
then for lower values of 7, reaching a = 0.86 in the limiting case 7 = 0.5. All 
our results for a, with 7 = 0.9, are thus still conservative estimates. 



3.1.2 Observational consequences 

We now switch to an observer's point of view. The observational data of a spe- 
cific cluster, like surface brightness and energy spectrum, do already contain 
the magnetic field effect! As argued before, fitting these data with an isother- 
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mal /3-model does not correspond to the real physical picture. However, the 
still large error bars in the present data will not allow to discriminate between 
an isothermal /3-fit and a fit with the profile pb- This is in particular true 
for the value of the central cluster density, which is often affected by errors 
of the order of 10% (see e.g. [Mohr et al.| ( |1999| )). The isothermal /5-profile 
PgiPpit) with the parameter Pfu, determined from observations, must thus be 
compared to p^: 

PB = p^m^ + hBir)] ^ PgiPFit). (35) 

As the central cluster region is often not very well resolved in X-ray data, Pfu 
is usually determined using the data from the outer parts. In these regions, 
the magnetic field influence is very small, i.e. /i^ ^ 0. In good approximation 
we get therefore: 

P ^ Pfu. (36) 



The isothermal /3-model will thus overestimate the real density[^ Pb{i^), be- 
cause hB{r) < 0. The error will be largest in the central region where the 
magnetic field pressure grows fast and becomes important. The central gas 
density can be overestimated by 10 — 20%, as it is seen in Fig. 1. We, there- 
fore, expect the observed SZ signal to be smaller by ~ 10% compared to the 
expected value as estimated from X-ray data, when the magnetic field influ- 
ence has not been taken into account. 

Following these arguments, we applied Eq.(ll) to a well documented non- 
cooling flow cluster: A119, which is characterized by the presence of three ex- 
tended radio galaxies. From different observations remarkably consistent val- 
ues are found ( |Mohr et al.| ( |1999| ), polag et al.| ( |2001| ) and references therein) 
and the existence of a cooling flow can be excluded for all practical purposes 
( [White et al.| , [19971 ). For an isothermal /3-model, the cluster parameter^ are: 



pp^^ = 0.56, 7 = 0.9, Bo = 7.5 pG, Tg = 5.92 ■ WK, Ug^ = 2.593 ■ lO'^ cm'^, 
Tc = 800 kpc, ri = 1550 kpc. From our discussion above, the quoted central 
value Ug^o should be overestimated. The other values are acceptable in good 
approximation, as explained. The results for the normalized profiles of Pb{t) 
compared to pg{l3, r) are shown in Fig. 2 and the calculated ratio a turns out 
to be a ~ 0.905. As Ug^ is overestimated, this is a conservative result and 
the real influence of the magnetic fields on the SZ-effect could be even a few 
percent larger. Reducing the central density by 10% gives a ratio of a = 0.895. 
If magnetic fields are neglected in A119, the SZ-effect is thus overestimated 
by ^ 10%. 



^ These arguments justify the use of the standard cluster values to find an estima- 
tion for pB{f) in section 3.1.1. 

The cluster 'geometrical' parameters are supposed not to be affected by the mag- 
netic field. 
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Fig. 2. The modified profile psir) compared to the /3-profile pg{f3,r) for A119. The 
mean molecular weight is p = 0.63. The profiles are normalized by the central gas 
density pg^. 

3.2 Cooling flow clusters 



Majumdar and Nath ( pOOOP already studied the effect of cooling flows on the 
SZ distortion and the possible cosmological implications, in particular on the 
Hubble constant. In this section we analyse how magnetic fields, in top of 
their result, can influence the SZ signal. First, we discuss the solution of the 
Eqs.(22) and (23) with magnetic fields, and second, we give the results of our 
SZ calculation for the magnetic cooling flow model described in section 2.2. 
We distinguish the cooling flow region from the region outside TcooI- The later 
one is analysed following the result obtained in section 2.1. 



Once the sonic radius is found from the Eqs.(28) and (29), integrating away 
from Ts requires then to find nonsingular expressions for the derivatives of the 
differential Eqs.(22) and (23). This is a non-trivial task with the additional 
magnetic field contribution. Furthermore, this procedure (shooting method) 
requires an iterative process to match the hydrostatic boundary conditions at 
Tcooi- If we want to separate the magnetic field influence, this gets even more 
complicated as we would have to match both pBi.'^cooi) and pg{rcooi) with the 
same cooling flow parameter m. We would thus have to find possible values 
for 9b{ts) and 9{rs) at the sonic radius - which differ by the magnetic field 
contribution in the Eqs.(22) and (23) - which must then lead to the required 
values for Pb{Jcooi) and Pgircooi) at the cooling radius. Since these values at 
Tcooi are derived by the procedure described in section 2.1, they are not in- 
dependent from each other neither, but related through the magnetic field 
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strength B{rcooi)- ^arazin and Whi"te| (|1987|) were already aware of how sensi- 
tive the integration of the Eqs.(22) and (23), in the absence of magnetic fields, 
can be. Not only three boundary conditions (p, 9, v) must be matched, but 
the possible values of them are additionally constrained by three physically 
motivated conditions like hydrostatic equilibrium, mass deposition rate and 
vanishing numerators in the Eqs.(22) and (23) for a regular remaining flow. 
Not surprisingly that the range of allowed boundary conditions to find a tran- 
sonic flow is extremely small. 

We stress that our goal is not to develop a sophisticated cooling flow model, 
but to get an estimate of how magnetic fields can modify the SZ-effect in a 
cooling flow cluster. We therefore avoid the time consuming and complicated 
procedure as outlined above, and start our integration with physically reason- 
able parameters from Tcooi towards r^. The transition at between subsonic 
and supersonic is always accompanied by shocks. Moreover, this most inner 
part will also be under the influence of a central galaxy. Since the interplay be- 
tween cooling flows, central galaxy, shocks and the SZ-effect is not clear, we do 
not attempt to find solutions to the Eqs.(22) and (23) inside the sonic radius 
Tg. We will, therefore, cut out this supersonic region for our SZ calculation. 
Converging cooling flows amplify the magnetic field strength enormously until 
the magnetic field pressure (-Pb) becomes comparable to the thermal gas pres- 
sure {Pg). Our numerical studies show that this occurs typically at ~ 10 fcpc. 
The ultimate fate of the magnetic field is still a matter of debate. Field line 
reconnection or convective motions resulting in buoyantly rising regions might 
be possible answers. Furthermore, [Christodoulou and Sarazin] (|1996|) pointed 
out, that more realistic magnetized models can in this region not be treated 
correctly in spherical symmetry. Based on the above arguments we decided 
to stop the integration of the Eqs.(22) and (23) at the pressure equipartition 
radius rB{> rg), where ^ ~ 1. 



We briefly discuss the input values for our results. The gravitating cluster mass 
A4{r) is fixed as described in section 2.2. We assume that the cooling flow gas 
makes no significant contribution to the mass density of the cluster. For the 
region outside rcooi, a standard isothermal /5-model with /3 = |, = 1.9 • 
10^ K and a central gas density o = 1-2 ■ 10~^ cm~^ is consistently adopted. 
We choose a magnetic field strength B{rcooi) ~ 1 /uG, which reproduces the 
observed values of the order of Bq k, 50 ixG for the cluster core, as found from 
the cooling flow model in the Eq.(20). Assuming again 7 = 0.9, these values 
determine then with the Eq.(ll) the density profiles pg{r) and psi'^) for the 
region outside Tcooi, and define the starting values Pg{rcooi) and pBiTcooi) for the 
cooling flow region. Once the mass deposition rate fn is fixed, the continuity 
Eq.(18) completely determines the initial velocities vb{tcooi) and v{rcooi) for 
the integration of the Eqs.(22) and (23). For these velocities at TcooI we must 
require that they are of the order of a few tens of km/ s, which corresponds to 
the turbulent velocity vt- From the difference between psircooi) and Pg{rcooi), 
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Fig. 3. The influence of magnetic fields on the coohng flow velocity \v{r)\ and 
the isothermal sound speed Cs{r) as a function of radius. The mass deposition 
rate is m = — 300 MQyr~^, r^ooi = lOOkpc and B{rcooi) = The initial 

conditions for the integration are: \/6{rcooi) = \'^h}(fcooi) = 1600A;m/s and 
\v{rcool) \ = \vB{rcool)\ = I5km/s. 

{pBircool) ~ Pgij-cooi)), it is obvlous that \vB{rcooi)\ ~ \v{rcooi)\- The initial 
sound speed squared, OBijcooi) = 0{rcooi), is directly related to Te = Tg. 
Finally, the so derived initial conditions and chosen parameters have to allow 
for the existence of a cooling flow: tcooi ~ tuubUe- From the above discussion it 
is clear that the integration with possible values from TcooI towards is also 
limited by constraints, which can be summarized as follows: 



M{r)=MB{r). (37) 

Tfl 1 

Pgircooi) = — 2 7 V' \'v{rcooi) \ ^ Vt, (38) 

^^ri^oi v{rcooi) 

Tfl 1 

PBircool) = — 2 7 T' \vB{rcooi)\ ^ Vt, (39) 

tcool = ~ ^Hubble ^ lO'^r. (40) 

2 pA 



These constraints have to be fulfilled with observationally reasonable values m 
and Bq. As a possible set of initial values and cluster parameters we have cho- 
sen: B{rcooi) = IpG, m = -300 Mo/yr, ^J6{r~i) = ^J9^{r~i) = 1600 fcm/s, 
\v{rcooi)\ = \vB{rcooi) \ = I5km/s. 

Fig. 3 shows how the cooling flow dynamics are influenced by magnetic fields. 
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Fig. 4. The influence of magnetic fields on the cooling flow gas number density. 
Parameters and initial conditions as adopted in Fig. 3. 

The integration for the profiles with magnetic fields is stopped at the pres- 
sure equipartition radius r^, whereas the profiles in absence of magnetic fields 
are stopped close to the sonic radius where the Mach number is M ^ 0.9. 
Fig. 4 and Fig. 5 give the corresponding gas number density and temperature 
profiles for the cooling flow region. Since the magnetic field strength at Vcooi is 
very weak, but the cluster temperature high, the difference between pg{rcooi) 
and pB{.fcooi) is completely negligible, as it results from Eq.(ll). Therefore, 
the initial velocities \v{rcooi)\ and \vb{j'cooi) \ can be set equal for all practical 
purposes. The influence of the magnetic fields on the SZ-effect is thus entirely 
determined by the cooling flow region. 

Similarly to Eq.(34), we define a ratio acF for the cooling flow {CF) region. 
For our parameters we find: 



where we cut out the most central part of the cluster and the lower inte- 
gration limit is taken to be the pressure equipartition radius = 10 /cpc. If 
the above ratio is completed with the SZ contribution from the isothermal 
region outside Tcoou (^s /rlo; Asol?^)'^?^, with piso an isothermal /3-model and 
Ti = 1000 kpc), one finds acF ~ 0.99. The presence of magnetic fields reduces 
thus only slightly the cooling flow correction to the SZ-effect. When calculating 
the ratio as between the magnetic field influenced profiles and the standard 



OiCF '■ = 




0.98, 



(41) 
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isothermal profiles, we find instead: 

._ ^Ib,cf _ IrTpBir)TB{f)df + Tj;i^^pUf)df _ ^^^^^ ^^2^ 
A/ Tg Jp^ pisoif)df 



compared to 1.02, which we get without taking into account the magnetic 
field in the cooling flow. This is not surprising because of the weakness of the 
initial magnetic field strength at Vcooi, which leads to practically identical ini- 
tial conditions. The magnetic field strength becomes only important towards 
the very center and there it then almost compensates the small cooling flow 
effect. For the chosen parameters, magnetic field and even cooling flow do not 
modify much the SZ-effect. 

Nevertheless, our numerical investigations showed that for a higher magnetic 
field strength, B{rcooi) ~ 6 /iG, and a higher temperature, ^ 5 • 10^ K, the 
ratio in Eq.(41) can reach acF = 0.94. Moreover, at variance with the above 
discussed example, a higher magnetic field strength at TcooI leads also to a size- 
able lower gas density, as compared to the case without magnetic field, in the 
region outside TcooI- This fact, as discussed in Section 3.1.1, modifies accord- 
ingly the contribution to the SZ-effect coming from the integration outside 
Tcooi- For such cooling flow parameters, jMajumdar and Nath| (|2000|) expect 
an overestimation of the standard SZ-effect by ~ 10%, due to the presence of 
a non-magnetized cooling flow. We instead find that this overestimation gets 
reduced to 5 — 8% when the (higher) magnetic field is taken into account. 
Finally we note that our results are conservative estimates since we always 
cut out the most central part of the cluster. 



4 Conclusions 



In a phenomenological approach we added the magnetic field pressure to the 
gas pressure and for clusters without cooling flows we derived in a perturbative 
procedure a new gas density profile. This can be related to a standard /3-profile 
and a function, which takes into account the radius dependent magnetic field 
strength, which is assumed to be correlated to the electron density. For rea- 
sonable cluster parameters we find that magnetic fields reduce the standard 
SZ signal by ~ 10%. Indeed, a reduction of up to ~ 15% seems possible. Our 
perturbative approach based on the equal mass assumption, 7VlB(r) = A1(r), 
turns out to be well justified by this order of magnitude correction: The corre- 
sponding decrease in the gas density is < 10 — 15% and, therefore, the change 
in the dark matter dominated total cluster mass profile is less than 1%. The 
simplifying assumption of an isothermal temperature is adequate, because the 
interplay between temperature and a variable magnetic field strength is not 
yet very clear ( polag et al.| , |2001| ) . 
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Fig. 5. The influence of magnetic fields on the cooling flow temperature. Parameters 
and initial conditions as adopted in Fig. 3. 

Furthermore, our considerations showed that the central cluster gas density 
is probably overestimated by 10 — 20% when fitted with a standard /?-model. 
Better data in future might reveal the need for a modified /^-profile as dis- 
cussed here. 

Other possible causes of deviations from the (standard) /3-model have been 
discussed in the literature in the context of the determination of the Hub- 
ble constant ( piirkinshaw et aT| ( |1991| ); [Inagaki et al.| ( |1995| ) and references 
therein). Generally, the finite extension of a cluster (already adopted in our 
calculation) lowers the SZ signal by 5 — 10% compared to the case of an 
infinite isothermal /3-model, and requires then a larger core radius Tc in the 
/3-model. Other departures from the standard /?-model include asphericity and 
small-scale clumping. The most extreme variation of geometry of the original 
spherical model is obtained if the unique axis of the prolate {Z > 1) or oblate 
{Z < 1) gas distribution is oriented along the line of sight. The SZ signal scales 
then with the factor Z, which is the ratio of the length of the unique axis to 
the major or minor axis, respectively. Typically, 0.5 < Z < 2, and asphericity 
introduces thus a variation of a factor ~ 2 in the core radius Vc- Whereas this 
effect will be averaged out in a large enough cluster sample, all cluster atmo- 
spheres will be clumpy to some degree. This might result in an overestimation 
of the standard SZ signal at the percentage level, but current simulations are 
limited in resolution and dumpiness might be more important. Whereas all 
these effects modify the parameters of the standard /?-model, the inclusion 
of magnetic fields requires an extended /5-profile, which takes explicitly into 
account the magnetic field strength. Unfortunately, in an integrated SZ mea- 
surement all the effects will sum up and they can hardly be separated. Precise 
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future X-ray data might then help to reveal the inner cluster structure. 



For the cooling flow clusters, following the treatment by |Soker and SarazTn] 
(|1990|) , we generalized the equations for the sound speed and infall velocity 
derived by [Mathews and Bregman] (|1978|) such as to include also a magnetic 



fleld. For typical initial magnetic fleld values at TcooU Pb/Pq ~ 10"'^, we flnd an 
equipartition radius ~ 10 kpc and we derive proflles which are in agreement 
with Soker and Sarazin. Our somewhat larger sonic radius ~ 5 kpc - though 
not irrealistic as it was found by pulkanen et al.| (|1989D , who derived sonic 



radii of ~ 10 — 20 kpc for plausible cluster and cooling flow parameters - 
results probably from different initial conditions. As the integration for the 
equations of sound speed and infall velocity is very delicate or even impossible 
for any arbitrary combination of parameters, our results can not immediately 
be generalized. Current discussions about the magnetic fleld model in the 
cooling flow region limit further our conclusions. Nevertheless, it turns out 
that the gas density in the cooling flow gets somewhat smaller in the presence 
of magnetic flelds as compared to the case without. This translates then into 
a weaker influence of the order of some percent of the magnetized cooling flow 
on the SZ-effect. In special cases magnetic flelds might almost compensate the 
cooling flow correction to the SZ-effect. 

Indeed, a precise calculation would require more sophisticated models which 
are adapted to the speciflc cluster parameters. For all these reasons, cooling 
flow clusters do not seem to be ideal targets for SZ observations. 
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